BMC Genomic Data

Chromatin profiling methods, such as ChIP-seq (chromatin immunoprecipitation followed by sequencing), are used to characterize the genomic localization of DNA-associated proteins. While ChIP and ChIP-seq have been used for decades, an orthogonal approach, CUT&Tag (Cleavage Under Targets and Tagmentation), is gaining popularity as an efficient and cost-effective alternative. Although previous comparative studies note discrepancies in signal-to-noise ratios and detection bias at certain genomic regions, many differences between ChIP-seq and CUT&Tag results remain largely unreconciled. Here, we systematically investigate the disagreeing signals captured by these two methods across well annotated genomic regions. We assess multiple histone mark profiles generated by different groups in two cell lines (K562 and MCF-7). Overall, our analysis indicates that compared to ChIP-seq, CUT&Tag may have limited sensitivity in low-GC environments and, as previously observed, increased signal in hyper-accessible chromatin. Notably, within GC-poor regions of active gene bodies and Polycomb-repressed domains, CUT&Tag exhibits a loss of H3K36me3 and H3K27me3 signal, respectively, where occupancy of these histone marks is otherwise expected. Further, active promoters with discrepant H3K4me3 and H3K27ac signal between the two assays differ systematically in GC content and chromatin accessibility. Promoters differentially enriched for CUT&Tag signal relative to ChIP-seq generally show higher, broader GC-content profiles and higher DNase-seq and ATAC-seq signal, while promoters enriched for ChIP-seq signal harbor the opposite characteristics. A local bias for high GC content and/or chromatin accessibility in CUT&Tag may also explain its differing signal patterns at nucleosome-depleted regions (NDRs) compared to ChIP-seq and MNase-seq promoter profiles. Altogether, our results highlight that studies focused on profiling the intensity and structure of histone modification occupancy can be sensitive to potential biases of CUT&Tag to GC content and chromatin accessibility. These characteristics should be accounted for when selecting a chromatin profiling approach, analyzing and interpreting data as well as drawing biological conclusions.

Human metabolism converts inorganic arsenic to the pentavalent methylated species MMA(V) and DMA(V), the forms most people excrete, and the forms long read as the end of a detoxification pathway. Whether a transporter sets how much of these metabolites reaches the genome has not been tested in a mammalian cell. We expressed human AQP3, AQP7, AQP9 or AQP10 in HEK293T and MRC5-SV40 cells and measured gamma-H2AX by flow cytometry across dose series of As(V), MMA(V) and DMA(V), pairing every aquaporin with a GFP-Tubulin control and an untransfected mock acquired in the same replicate. As(V) was inactive in HEK293T cells and only weakly active in MRC5-SV40 cells to 20 micromolar, and both methylated species damaged DNA only in the millimolar range, DMA(V) being the more potent of the two in both cell lines. Against that weak baseline, AQP9 and AQP10 raised DMA(V)-induced gamma-H2AX in HEK293T cells by roughly 17 percentage points over the matched control, more than doubling the damage the same exposure produced in control cells, whereas AQP3 and AQP7 changed it not at all. AQP9 alone remained active with MMA(V). The ranking held in MRC5-SV40 fibroblasts at one-sixth the size, and within single wells the damage rose with the amount of AQP9 a cell carried while the control was flat. Aquaglyceroporins therefore discriminate among arsenic species, and AQP9 and AQP10 turn a weakly genotoxic metabolite into a substantially more genotoxic one.

Chronic exposure to inorganic arsenic affects more than 200 million people, and disease outcome varies widely among individuals with comparable intake. Explanations have centered on arsenic metabolism, which acts on arsenite that is already inside the cell. Uptake sits upstream of metabolism and sets how much arsenite reaches the genome, yet whether it acts as a gene-by-environment modifier of genotoxicity has been asserted far more often than it has been measured. We expressed human aquaporin-3 (AQP3), its point mutants and the paralogues AQP7, AQP9 and AQP10 in three human cell lines, and measured DNA damage as {gamma}H2AX by flow cytometry against a damage threshold re-derived from mock-transfected cells in every experiment. In lung fibroblasts, AQP3 raised the {gamma}H2AX-high fraction 2.36-fold at 5 M arsenite. Because a gene-by-environment interaction widens the genotype gap with dose while an independent genotoxin would not, the combination exceeds what the two factors would produce if they simply added by 6.36 percentage points, which is an interaction. The response required the AQP3 sequence, was suppressed by N-acetylcysteine in proportion to arsenite dose, was absent in RKO cells, and was shared with AQP9 and AQP10, which created an arsenite dose-response in a cell line that had none of its own. Low-input RNA sequencing of 135 ten-cell pools located the effect inside the cell: at an identical 1 M applied dose and within a single batch, AQP3 raised HMOX1 (induced when trivalent arsenic modifies KEAP1 thiols, and so a report of intracellular rather than applied exposure) 3.3-fold, a shift equivalent to 40% of the 1-to-10 M interval on the measured dose-response, while AQP3 without arsenite was indistinguishable from untreated. Duplex sequencing detected no consistent change in somatic mutation frequency from arsenite alone at any dose to 10 M, in an assay that resolved a 26-fold response to N-ethyl-N-nitrosourea; the arsenite-induced mutation frequency at the EGFR locus nonetheless differed 1.61-fold by AQP3 genotype. The interaction is therefore not unique to AQP3, and aquaglyceroporin complement belongs alongside metabolism among the host factors that set arsenic susceptibility.

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants with the potential to disrupt metabolic processes and contribute to carcinogenesis; however, molecular evidence linking alternative and emerging PFAS to ovarian cancer (OC) remains limited. We conducted a biomarker-based case-control study integrating exposureomics with multi-omics profiling. Plasma samples from 128 patients with OC and 128 controls were analyzed for 26 PFAS, untargeted metabolomics, lipidomics, and clinical phenotypes, while urine samples from an additional 100 patients and 100 controls underwent untargeted metabolomic and clinical profiling. Single-compound, burden-based, and mixture models were integrated with PFAS-omics association scans, cross-matrix metabolite prioritization, observational exposure-omics-clinical association-chain analyses, and repeated internal cross-validation. Higher PFAS burdens were consistently associated with OC case status, and DONA, GenX, PFNA, PFPeS, PFUdA, and PFTrDA constituted an OC-prioritized six-PFAS mixture. DONA and GenX exhibited the broadest lipid-related bridging signals, particularly across phospholipids, lysophospholipids, sphingolipids, neutral lipids, and acylcarnitines. Integrated plasma and urine analyses further identified five shared metabolites-elaidic acid, hippuric acid, betaine, corticosterone, and taurocholic acid-that provided incremental and complementary discrimination beyond reference models based on clinical variables and tumor markers in internal cross-validation. Collectively, these findings demonstrate that both alternative and emerging PFAS are associated with a distinct OC-related molecular phenotype involving coordinated lipidomic, metabolomic, and hematologic alterations, while highlighting DONA and GenX as the most prominent exposure-related signals.

Fetal development is highly sensitive to environmental influences, including maternal smoking during pregnancy (MSDP) and air pollution. The placenta may play a central role in mediating these effects through epigenetic mechanisms, including small non-coding RNAs (sncRNAs). While most studies have focused on microRNAs (miRNAs), other sncRNA classes remain poorly characterized. Here, we characterized the placental sncRNA landscape and investigated its associations with prenatal exposure to MSDP and air pollution. The sncRNA expression profile of 390 chorionic villi samples from the Barcelona Life Study Cohort (BiSC) was measured using small RNA sequencing and the seqBuster/seqCluster pipeline. MSDP was assessed via questionnaire, and prenatal exposure to black carbon (BC), particulate matter [≤]2.5 m (PM2.5), and nitrogen dioxide (NO2) was estimated using hybrid dispersion-land-use regression (LUR) spatiotemporal models. Associations were evaluated using linear regressions adjusted for relevant covariates. We identified 979 placental sncRNA clusters, assigned primarily to miRNAs (33.2%), piRNAs (26.4%), snoRNAs (12.2%), tRFs (7.8%), and others (20.4%). MSDP was associated with the downregulation of miR-421 and the upregulation of multiple snoRNAs from the imprinted DLK1-DIO3 locus (SNORD113-4, SNORD114-2, SNORD114-3, and SNORD114-17), as well as the piR-26394. No sncRNAs remained significantly associated with air pollution after multiple-testing correction, although several miRNAs showed suggestive associations. Overall, this study provides a comprehensive characterization of placental sncRNAs and identifies previously unreported smoking-associated sncRNAs within the imprinted DLK1-DIO3 locus, expanding our understanding of placental epigenetic responses to prenatal environmental exposures.

Vertical transmission of SARS-CoV-2 is controversial. Owing to the low expression of ACE2 along with TMPRSS2 in placenta, role of secondary receptors becomes significant. In our study we predicted twelve interactors of spike protein expressed in placenta through interactome analysis followed by docking and molecular dynamic simulation. Three interactors of SARS-CoV-2 spike protein were found to be CD147, EGFR and LDLR. Expression at transcriptomic and translational level by real-time PCR and immunohistochemistry respectively was observed in placenta samples obtained from women with COVID exposure risk and compared to control. Elevated protein levels of CD147 (OR=14.667), ACE2 (OR=8.25) and LDLR (OR=1.9) was associated with higher exposure risk. However, the variation in mRNA expression was not statistically significant. Computational study showed that CD147/basigin, EGFR and LDLR stably interact with spike protein outside its primary receptor binding domain. In-vitro binding assay confirmed the physical interaction between spike protein and the placental receptors. In-silico site-directed mutagenesis further identified key potential binding hotspots of spike protein- Tyr28 and Trp64 for CD147, Phe643, Val642, Thr632 for LDLR and Leu8, Phe59, Phe32, Tyr28, Trp64, Trwp152, His245, and Val213 for EGFR. Finally, concluded that CD147, EGFR and LDLR might act as secondary non-canonical cofactors for SARS-CoV-2 spike protein in placenta, thereby increasing its susceptibility and tissue tropism. The study overall gives a better insight into potential placental suscpetibility providing targeted therapeutic opportunities.

CBX1-related syndrome is characterized by developmental delay, hypotonia, autistic features, and mild dysmorphic features. This syndrome is caused by heterozygous missense variants in CBX1, which encodes heterochromatin protein 1 beta (HP1{beta}). These variants are situated within the chromodomain of HP1{beta}, a critical region that mediates the interaction between HP1{beta} and trimethylated H3K9. While a prior study identified three individuals with pathogenic CBX1 variants, the variant and clinical spectrum of CBX1-related syndrome remains largely unknown. Here we report three patients with novel de novo heterozygous variants in CBX1. Patients with CBX1 p.(Leu40Pro) and p.(Cys60Arg) manifested with global developmental delay, while the patient with CBX1 p.(Cys60Gly) primarily exhibited speech delay without concurrent motor developmental delay. All variants reside within chromodomain of the HP1{beta} protein. Interestingly, in vitro cellular assays revealed that CBX1 p.(Cys60Gly) variant showed an intermediate cellular phenotype between the wild type and the CBX1 p.(Leu40Pro) and p.(Cys60Arg) variants, suggesting a possible genotype-phenotype correlation. This report provides further evidence that missense variants within the chromodomain of HP1{beta} cause a syndromic neurodevelopmental disorder, CBX1-related syndrome.

In eukaryotic cells, the positive transcription elongation factor b (P-TEFb) plays a critical role in the transition of RNAPII from the paused state to actively transcribing mode. Although three different cyclin (Cyc) T (CycT1, T2a, and T2b) interact with CDK9 to form functional P-TEFb, a majority of P-TEFb complexes are comprised of the CycT1 and CDK9. To analyze the role of the CycT1 subunit of P-TEFb on cellular transcription, we established HEK293T cells lacking CycT1 proteins (CycT1-KO cells). While no apparent growth defects were observed with CYcT1-KO cells, the level of CDK9 was decreased and no compensatory over-expression of CycT2 was observed, suggesting that the reduced level of functional P-TEFb complexes is sufficient for the normal cell growth. Reporter gene assays indicate that NF-kappaB-dependent transcription induced by PMA was impaired in CycT1-KO cells although the function of NF-kappaB (P65) per se was not affected. On the other hand, AP1-dependent transactivation induced by PMA was unaffected in CycT1-KO cells. Moreover, transcription stimulated by JQ1, a strong P-TEFb inducer, was reduced in CycT1-KO cells. Data of biochemical analysis indicate that P-TEFb with CycT1 (P-TEFb (CDK9:CycT1)) was efficiently released from 7SKsnRNP by JQ1 while P-TEFb (CDK9:CycT2) was not affected by JQ1. Transcriptome analysis indicated that the lack of CycT1 had a minor effect on the steady-state transcription. However, expression of JQ1-dependent genes was severely reduced in CycT1-KO cells. Interestingly, referring to gene expression profiles of 33 different types of cancer revealed that 11 out of 13 top JQ1-dependent genes were also highly upregulated in Pancreatic Adenocarcinoma. From these results, we conclude that although many genes can be regulated by both P-TEFb (CDK9:CycT1) and P-TEFb (CDK9:CycT2), CycT1 plays a critical role in regulating highly inducible genes, which are also aberrantly regulated in a particular type of cancer.

Heavy metal contamination of aquatic environments poses persistent risks to ecosystems and human health. Aquatic macrophytes such as Pistia stratiotes are increasingly used for phytoremediation due to their capacity to accumulate contaminants without displaying symptoms of abiotic stress. However, the molecular mechanisms underlying their metal tolerance remain poorly understood. Metallothioneins (MTs) are cysteine-rich proteins involved in metal binding and detoxification in plants. In this study, genome-wide analysis was performed to identify and characterise candidate MT genes in P. stratiotes and to assess their transcriptional response to copper exposure. Six novel MT genes (PsMT2f, PsMT3a-d, and PsMT4) were identified and classified into MT2, MT3, and MT4 subfamilies based on conserved cysteine motifs and domain architecture. Quantitative real-time PCR showed differential and tissue-specific expression of MT genes in roots and leaves following exposure to copper for up to 4 h. These findings provide the first comprehensive characterisation of the MT gene family in P. stratiotes and demonstrate their potential role in copper tolerance, supporting the application of this species in the phytoremediation of metal-contaminated aquatic systems.

Individuals carrying a G4C2 repeat expansion in C9orf72 can develop amyotrophic lateral sclerosis (ALS) and/or frontotemporal dementia (FTD), at various and unpredictable ages. The contribution of the expansion to this extreme clinical variability remained a matter of debate. In this work, we validated the use of optical genome mapping to efficiently size this expansion from blood and applied this method to a cohort of 54 repeat expansion carriers including patients with FTD (n = 23), ALS (n = 22) and presymptomatic individuals (n = 9). Despite the high variability in expansion profiles, distinct groups can be inferred. Although consensus repeat profiles were similar between FTD and ALS patients, relationships between repeat number, indexes of repeat instability, and the age at sampling differed between clinical conditions. An age at sampling-dependent increase of repeat instability in blood was observed in presymptomatic carriers and FTD patients (p = 0.01 and p = 0.04 respectively), in contrast to patients with ALS (p = 0.9). The proportion of shorter repeats ([500-1500 repeats]) increased with the age at sampling in ALS patients but decreased in presymptomatic and FTD cohorts, while the proportion of larger repeats (>2000 G4C2), in particular those exceeding 4000 repeats, increased along with the age at sampling in presymptomatic carriers and FTD cases (p <1.10-e4 and p = 0.006 respectively) conversely to patients with ALS (p = 0.98). In conclusion, this work strongly supports different dynamics of repeat instability of the C9orf72 pathogenic expansion in blood linked to the clinical evolution. These findings provide new insights into the relationship between expansion profiles and clinical variability.

Background: Across phenotypic and genotypic studies, executive functioning (EF) has been shown to be associated with externalizing behaviors, including attention-deficit/hyperactivity disorder (ADHD), aggression, and substance use, and is a proposed underlying mechanism of these behaviors. However, there is limited research on the shared biological mechanisms underpinning EF and externalizing behaviors, including their overlap at the level of genomic regions. Methods: The current study used genome-wide association study summary statistics of EF, ADHD, aggressive behavior, alcohol use, and alcohol use disorder to examine whether the heritability of these traits or the co-heritability between them is enriched for genes expressed in particular tissues or stratified by other annotation types. Additionally, the study examined the genomic regions contributing to the covariation between EF and the externalizing traits to identify genetic variants and genes that may be underlying their overlap. Results: First, across all phenotypes, there was evidence for enrichment of genetic signals in DNase hypersensitive sites and H3K4me1 regions, and significant enrichment among genes expressed in the anterior cingulate cortex and frontal cortex. Second, when combining results across regional methods, there was evidence for associations between EF and each of the externalizing behaviors in regions on chromosomes 1, 11, and 15. Conclusions: These findings highlight biological pathways that may link the relations between EF and externalizing behaviors and provides additional evidence of EF as an underlying mechanism of externalizing behavior.

Abstract Polycystic ovarian syndrome and Type 2 diabetes mellitus are complex multivariate diseases that share well-established clinical pathophysiology, yet the molecular basis of their overlap is unknown. Prior bioinformatics studies have identified common differentially expressed genes between these two conditions. However, they have not accounted for the directionality of the gene expression changes observed, potentially obscuring any biological distinctions they may have. This study applied a direction-aware approach to classify molecular signatures between PCOS and T2DM as concordant or discordant. Transcriptomic analysis of GEO datasets GSE138518 (ovarian granulosa tissue, PCOS) and GSE25724 (pancreatic islet tissue, T2DM) identified 225 and 1,302 DEGs, respectively. Venn diagram analysis showed that only three genes (SLC6A8, RGS4, and SORL1) were shared between PCOS and T2DM from the total of 1,527 genes, and all 3 genes were regulated in opposing directions. Disease gene retrieval from the Comparative Toxicogenomics Database, Online Mendelian Inheritance in Man database, and GeneCards showed 214 shared disease-associated genes, with 311 genes unique to PCOS and 496 genes unique to T2DM. Protein-protein interaction construction using STRING (version 12.0) identified 70 interacting nodes. CytoHubba analysis across 6 scoring methods identified 18 high-confidence hub genes, including INS, BCL2, MTOR, LEP, MFN2, and PIK3CD. Functional enrichment analysis identified biological processes including immune activation, apoptosis, and insulin signaling, which were confirmed as the main pathways in KEGG pathway analysis. These findings show that while PCOS and T2DM share limited DEGs, they share common pathogenic pathways.

Background: Preeclampsia (PE) is a hypertensive pregnancy disorder with contributions from both maternal and fetal genetic factors, but how both genomes regulate placental transcription remains incompletely understood. Results: We analyzed placental 5'-end RNA sequencing data from 59 pregnancies in the FINNPEC cohort and performed separate maternal and fetal cis-eQTL analyses using gene-level expression and promoter-level expression defined as Transcript Far 5' Ends (TFE). Genotype-by-PE interaction models were used to identify genetic effects that differ according to PE status. No gene-level interaction eQTLs remained significant after multiple-testing correction, although suggestive signals included LILRB2 in the maternal and LPCAT3 in the fetal analyses. In contrast, we identified two significant syntenic promoter-level interaction eQTLs: rs12504642-TFE61820 at the FHDC1 locus in the maternal analysis and rs1557150-TFE9032 at the PPIAP31 locus in the fetal analysis. Conclusion: These findings suggest that PE-dependent genetic interaction effects may be detectable at specific transcript promoter regions even when they are not captured at annotated gene levels.

Background: Long non protein coding (lnc)RNAs are ubiquitous transcripts known to be major regulators of cellular processes in eukaryotes. In diatoms, specific lncRNA expression is tightly correlated to nutrient fluctuations, which suggests function, but their biological relevance remains to be established. In order to fill this knowledge gap, we sought to design a protocol to generate diatom lncRNA knockout mutants by CRISPR-Cas9 enabling subsequent functional studies to shed light into their mechanistic roles. Since lncRNAs do not contain open reading frames that can be targeted by point mutations for the generation of functional knockouts, we designed a multiplexed CRISPR-Cas9 plasmid harboring two independent gRNA for the specific ablation of lncRNA loci in P. tricornutum. As proof of concept, we used our plasmid-based system to target three phosphate depletion responsive lncRNAs. Results: Using the golden gate based uLoop system, we assembled a single plasmid containing all the necessary components to generate CRISPR-Cas9 large genomic ablations of lncRNAs in P. tricornutum using bacterial conjugation. Three intergenic lncRNAs (Linc91, Linc77 and Linc24) responsive to phosphate fluctuations were selected and independently targeted. Around 32% of the exconjugants presented a monoallelic ablation of the targeted loci and 11% presented a biallelic ablation on the first PCR-screening round. Mutated lines, harboring biallelic lncRNAs ablations for each lncRNA targeted, were successfully generated for each of the targeted lncRNAs. Additionally, a first line of phenotypic analysis revealed significant growth rate differences for one of the mutants ({Delta}Lnc24) under both full and phosphate depleted conditions, suggesting its involvement in cell division regulation under nutrient fluctuations. Conclusions: This work provides a simple and efficient method for the generation of lncRNA knock-out mutants in P. tricornutum. Furthermore, this approach can be easily adapted to other microalgae amenable to transformation, enabling the enhancement of our understanding of these important cellular regulatory transcripts in the marine environment.

Background: More than five years after the COVID-19 outbreak, the infection by the SARS-CoV-2 virus should still be considered as a global threat with highly heterogeneous clinical outcomes from asymptomatic to fatal depending on existing genetic and phenotypic inter-individual and inter-population variability. Our objective is to identify risk loci associated with COVID-19 in an admixed Latin American population. Methods: We performed a genome-wide association study in the framework of the Chilean C19-GenoNet team, analyzing a cohort of 3,417 individuals recruited from the northernmost to the southernmost regions of the country, who presented variations in the severity of their COVID-19 clinical manifestations. Results: Our investigation replicated the extensively reported association of several variants at the 3p21.31 locus with COVID-19 severity and found a suggestive association for a variant within the FOXP4 gene at 6p21.1 with COVID-19 hospitalization. Both associations appear to be shared with other populations. Importantly, our study identified three previously unreported loci associated with COVID-19 severity or hospitalization at 7q21.11, 16p13.12 and 6p25.3, and one additional locus at 4p13 associated with SARS-CoV-2 positivity, a proxy for infection susceptibility. Conclusions: these new risk loci may represent population-specific effects in Chilean and, potentially, other Latin American admixed groups. Our study highlights the importance of expanding efforts to study underrepresented and genetically different populations, which may provide new clues on the genetic basis of human diseases and may help to anticipate global responses to emerging viral infections and future pandemics.

Minor introns are rare but evolutionarily conserved introns in the human genome. Defects in minor spliceosome components have been reported to cause several human diseases. Mutations in ZRSR2 are frequently observed in myelodysplastic syndrome (MDS) patients. However, the precise molecular mechanism by which these mutations cause a loss of ZRSR2 function in minor intron splicing remain elusive. In this study, we characterized the impact of MDS-causing missense mutations in ZRSR2. Most of the mutant ZRSR2 proteins we tested showed different subnuclear localization pattern from that of wild-type ZRSR2 protein. Co-immunoprecipitation (co-IP) analysis revealed that UHM missense mutations disrupted the interaction between ZRSR2 and SF3B complex. We also demonstrated that ZRSR2 preferentially immunoprecipitated minor spliceosome components over major spliceosome components and it directly interacts with SF3B1. Furthermore, co-IP analysis utilizing a variety of ZRSR2 mutants revealed that mutations in the UHM and Zn2 domains significantly attenuated the interaction with SF3B complex. Finally, UHM and Zn2 mutations markedly disrupted the supporting activity of ZRSR2 for minor intron splicing. Our results indicate that disruption of the physical interaction with SF3B complex is one of the hallmarks of MDS-causing mutations of ZRSR2, which results in aberrant splicing of minor introns.

Little is currently known about the regulation of PL fimbriae produced by some pathogenic E. coli strains. We identified PlfB, a transcriptional regulator of the PapB family, as a negative regulator of PL fimbriae (Plf). Deletion of plfB significantly increased plfA transcription and Plf production. PlfB bound specifically to the plfA promoter region. Using a PplfA-lux reporter system, the deletion of plfB caused sustained and increased expression across serial passages. Global regulatory genes such as lrp, argR, and hns were also implicated in regulation of plf expression. Loss of lrp or argR decreased plf expression, suggesting these regulators positively regulate PL fimbriae. By contrast, loss of hns increased levels of plf expression, supporting its role as a negative regulator. Interestingly, loss of plfB significantly increased bacterial colonization in the lung of 6-day-old turkey poults. Further, complementation of plfB restored plf fimbrial expression and colonization of turkey lungs to wild-type levels. This further supports an important role for PL fimbriae for colonization of turkey lungs, since loss of plfB which resulted in increased plf gene expression in turkey lungs also correlated with higher bacterial numbers in the lungs. Taken together, these results indicate that PlfB acts as a negative regulator to reduce levels of PL fimbrial expression, under some conditions, but that it is not required for plf expression. This contrasts with the role the PapB regulator, which has been shown to play an essential role in expression of P fimbriae in uropathogenic E. coli strains.

Abstract Background. Tuberous Sclerosis Complex (TSC) is caused by inactivating mutations in either the TSC1 or TSC2 genes, leading to activation of the mammalian target of rapamycin complex 1 (mTORC1) and unhindered cell growth and proliferation. The epithelium of TSC renal cysts in both mice and humans is composed of proliferating A-intercalated (A-IC) cells. The exact molecular mechanism of kidney cystogenesis in TSC remains speculative. Hypothesis.Superfluous cell proliferation driven by mTORC1 activation increases metabolic demand, causing oxidative stress and excess reactive oxygen species. If unchecked, this overwhelms antioxidant defenses and causes cell death. In TSC kidney cystogenesis, Nuclear Factor Erythroid 2-Related Factor 2 (NFE2L2), also known as NRF2, serves as the "master regulator" of antioxidant and anti-inflammatory responses, enabling cells to survive and proliferate by clearing the toxic environment and supplying nutrients and fuel. Results. RNA-seq and proteomics, along with western blot analysis, showed robust downregulation of Fumarate Hydrase 1 (FH1) and upregulation of NRF2, STAT3, and HIF1 in kidneys from TSC mice with moderate or heavy cyst burden. Confocal microscopy and immunohistochemical staining on kidney sections, and/or western blot studies on nuclear and cytoplasmic fractions, showed nuclear localization of NRF2, STAT3, and HIF1. In cyst-lining cells in TSC mouse models, FH1 downregulation was associated with inactivating succination of KEAP1 in immunoprecipitation experiments, promoting NRF2 nuclear localization in A-IC cells lining the cysts. NRF2 nuclear localization was associated with ectopic induction of the glutamine transporter SLC38A3 (SNAT3) on the basolateral membrane and activation of the NH3/NH4+ transporters RHCG and RHBG in A-IC cells lining the cysts. Twenty-four h urine NH3/NH4+ excretion rates increased significantly in Tsc1 KO vs. WT mice. The activation of NRF2, STAT3, and HIF1 can drive metabolic reprogramming and activate survival genes in proliferating cells. Together with SLC38A3 induction and upregulation of other glutamine and NH3/NH4+ transporters, these factors activate glutaminolysis and aerobic glycolysis, supplying nutrients to proliferating cystic epithelial cells and supporting cyst expansion in TSC. Consistent with this central role for glutaminolysis in kidney cystic epithelium and TSC cystogenesis, we find a significant reduction in kidney cyst burden in Tsc1 KO mice on a glutamine-free diet. Conclusions. NRF2 plays a critical role in antioxidant defense. Along with STAT3 and HIF1, NRF2 is a key player in metabolic reprogramming through glutaminolysis, which supplies nutrients to proliferating cystic epithelial cells and supports cyst expansion in TSC. These findings suggest that inhibiting or inactivating NRF2, alone or in combination with HIF1 or STAT3, may represent a potential treatment strategy for kidney lesions in TSC.

High throughput sequencing necessitates robust computational tools to contextualize complex transcriptomic data. In this study, we analyzed RNA Seq data from peripheral blood of bacterial pneumonia patients and healthy controls (GEO GSE161731) to investigate systemic molecular responses to infection. Differential gene expression analysis identified significantly upregulated genes in pneumonia. We used GANGO for Gene Ontology (GO) term enrichment and BioFunctional for hierarchical exploration, providing an initial overview of affected biological processes. To deepen interpretation, we leveraged CURIE (www.curieai.tech ), an AI plat-form that cross-references GO terms with scientific literature. CURIE generates narrative explanations, link-ing molecular mechanisms, experimental evidence, and publications to each biological process, moving be-yond simple enrichment. Our combined analysis revealed a coordinated upregulation of pathways involved in DNA repair and replication, transcriptional activation, RNA processing, translation, ribosome biogenesis, mitochondrial respiration, and cofactor binding. By integrating these findings within a broad, literature-supported framework, we highlighted how multiple cellular compartments and processes collaborate to mount an effective response to infection-induced stress. This study underscores the potential of combining classical statistical enrichment with AI based interpretation for a holistic understanding of host-pathogen in-teractions at the transcriptomic level.

Two rounds of whole-genome duplication (2R-WGD) in the vertebrate lineage are thought to have driven gene diversification. Mammals have four genes encoding five nebulin superfamily proteins (LASP1, LASP2/NEBL, NEB, and NRAP) hypothesised to have arisen from invertebrate LASP via WGD. Here, we examined the evolutionary fate of the nebulin superfamily after an additional WGD, teleost-specific genome duplication (TGD). Genome and transcriptome analysis of zebrafish (Danio rerio) and medaka (Oryzias latipes), of which the complete genomes were recently determined, showed only lasp family genes retained duplicate copies. Zebrafish retained two LASP2/NEBL-derived ohnologs; one (LASP2L) had been misidentified as a nebulette gene in the Zebrafish Information Network (ZFIN) despite lacking nebulette-specific exons. Another (LASP2/NEBL) is newly identified in this study because it lies in a newly determined region of the latest genome assembly. Contrastingly, medaka retained two LASP1 ohnologs and one LASP2/NEBL. The duplicated lasp family proteins showed overlapping and specific tissue expression patterns. We confirmed all three zebrafish lasp proteins bound actin and co-localized with actin-rich structures when ectopically expressed in cultured fibroblasts. This supports that TGD-derived lasp duplicates underwent lineage-specific retention and expression divergence in teleosts whilst conserving core actin-binding properties.